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Abstract 



Free-space propagation can be described as a shearing of the Wigner distribution func- 
\£) • tion in the spatial coordinate; this shearing is linear in paraxial approximation but assumes 

a more complex shape for wide-angle propagation. Integration in the frequency domain al- 



lows the determination of near-field diffraction, leading to the well known Fresnel diffrac- 
tion when small angles are considered and allowing exact prediction of wide-angle diffrac- 
tion. The authors use this technique to demonstrate evanescent wave formation and diffrac- 
tion elimination for very small apertures. 

The Wigner distribution function (WDF) provides a convenient way to describe an optical 
signal in space and spatial frequency ID E| . The propagation of an optical signal through first- 
order optical systems is well described by the WDF transformations 001, which allows the 
reconstruction of the propagated signal. On the other hand, some authors have linked Fresnel 
diffraction and the fractional Fourier transform (FRFT) 10 0; both of these papers associate 
free-space propagation to a rotation of the WDF and the corresponding FRFT accompanied by 
a quadratic phase factor. 



In this paper we show that free-space propagation is always associated with shearing of the 
WDF. This can be used to evaluate the near-field diffraction. In the paraxial approximation our 
results duplicate the well known Fresnel diffraction, with the advantage that we do not need to 
resort to the Cornu integrals. The same procedure can be extended to wide angles, where other 
phenomena are apparent, namely the presence of evanescent waves. 

The Wigner distribution function (WDF) of a scalar, time harmonic, and coherent field 
distribution <£>(q, z) can be defined at a z = const, plane in terms of either the field distribution 
or its Fourier transform ^(p) = / y?(q) exp(— ikq T p)dq J2j|3|: 

W(q, p) = J ip ^q + ^jf* (q - ex P (-^q' T p) d q (1) 

= / ^ (p + 2") ^* ( P ~ 2") 6XP ( ikqTp> ) dp ' ' (2) 

where q is the position vector, p the conjugate momentum, k = 2n/X and * indicates complex 
conjugate. 

In the paraxial approximation, propagation in a homogeneous medium of refractive index n 
produces a perfect mapping of the WDF according to the relation 

W(q,p,*)=W(q--p,p,0). (3) 

n 

After the WDF has been propagated over a distance, the field distribution can be recovered by 
O 

(p{q, z)<p*{0, Z ) = 2^J W W 2 > P' z ) exp(iqp)c/p. (4) 

The field intensity distribution can also be found by 

4tt 2 r 

Mq,z)\ 2 = —j W(q,p,z)dp. (5) 

Eqs. © and © are all that is needed for the evaluation of Fresnel diffraction fields. Con- 
sider the diffraction pattern for a rectangular aperture in one dimension. The WDF of the 



2 



aperture is given by 

\q\ > 1/2 

W(q, p) = { 2 sin[kp{l - 2q)}/kp 0<q<l/2 , 
2 sin[kp(l + 2q)]/kp -l/2<q<0 
with / being the aperture width. After propagation and integration in p we obtain 



fc2 M g '* )|2 = ^-sm[kp(l-2zp/n-2q)]dp 

8ll 2 Jn{2q~l)/(2z) Kp 

rn(2q+l)/(2z) \ 

+ / — smlkpCl + 2zp/n + 2q)]dp. 

Jnq/z kp 



(6) 



(7) 



Fig-dshows a typical diffraction pattern obtained in this way. 

For wide angles the WDF can no longer be mapped from the initial to the final planes. It 
must be born in mind that the components of the conjugate momentum p correspond to the 
direction cosines of the propagating direction multiplied by the refractive index; the appropriate 
transformation is given by Q 



W(q,p,z) 



W(q — zp/ yn 2 




|p| 2 ,P,0) |p| < n 
otherwise 



(8) 



Eq. ® shows that only spatial frequencies corresponding to momenta such that |p| < n 
can be propagated 0. In fact |p|/n = sin a, with a the angle the ray makes with the z axis. 
It is then obvious that the higher frequencies would correspond to values of | sin a\ > 1; these 
frequencies don't propagate and originate evanescent waves instead, Fig. [2j The net effect on 
near-field diffraction is that the high-frequency detail near the aperture is quickly reduced. 

The field intensity can now be evaluated by the expression 



k 2 \<p(q,z)\' 
4tt 2 



/ — sin <. kp 

Jpi kp I 

f P2 1 • f, 
+ / — sm < kp 

J vn kl) L 



I — 2 ( q — zp/ Jn 2 — p 



dp 



I + 2 ( q — zp/ J n 



V 



dp, 



(9) 
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with 



Po 




Pi 



n(2q - /)/^(2g-/) 2 + 4z 2 , 



n(2q + l)/^(2q + If + Az 2 . 



(10) 



Fig. |3] shows the near- field diffraction pattern when the aperture is exactly one wavelength 
wide. The situation is such that all the high order maxima of the WDF appear at values of 
\p\ > 1 and are evanescent, resulting in a field pattern with one small minimum immediately 
after the aperture, after which the beam takes a quasi-gaussian shape, without further minima. 
A sub-wavelength resolution would be possible about half a wavelength in front of the aperture, 
where the field distribution shows a very sharp peak. 

We return now to the subject of FRFT to analyze the connection between this and Fres- 
nel diffraction. Being associated with a rotation of the WDF, the FRFT can only represent 
diffraction if associated with a quadratic phase factor that effectively transforms the rotation 
in a shearing operation 0. The implementation of the FRFT needs a combination of paraxial 
free-space propagation, a thin lens and a magnification telescope JUIS). 

J. B. Almeida wishes to acknowledge the fruitful discussions with R Andres, W. Furlan and 
G. Saavedra at the University of Valencia. 
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Figure 2: Propagation of the WDF in wide angle condition (k = 10 6 m 1 , horizontal scale in 
jum). a) Original distribution, b) after propagation over 3 fxm. 




Figure 3: Near- field diffraction pattern when the aperture width is exactly one wavelength; 
(k = 10 6 m _1 , both scales in /im). 
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